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Micromagnetic simulation results on Co-based recording media are presented which examine the impact of a modified near- 
neighbor exchange interaction between grains of the form / z (M[ Z Mj Z ), reflecting the hexagonal crystal symmetry. Both out-of-plane and 
in-plane M-H loops are calculated, with an emphasis on a model fit to data reported by Wang et al. [IEEE Trans. Magn. vol. 43, 682 
(2007)] on exchange coupled composite perpendicular media. The principle effect of J 1 is to increase the coercivity and slope of both 
hard and soft layers. Improved agreement with experimental data for the in-plane loops is achieved by assuming a substantial value 
for f\ Possibilities for measurement of J z through spin-wave excitations are discussed. Thermal fluctuation effects are also examined 
through simulations of the magnetization vs temperature. 

Index Terms — Micromagnetic, perpendicular media, hysteresis, composite. 



I. Introduction 

The amorphous nature of sputtered cobalt-based 
granular thin films for recording media offers a number of 
challenges to model key features relevant to increasing areal 
density and thermal stability[ 1] . The fundamental ingredient 
to modern micromagnetics is the energy expressed as a 
functional of the magnetic moment vectors of interacting 
grains, E[MJ. An important aspect of the energy is that it 
should reflect the symmetry of the granular crystallites, 
typically hep for Co-based films. It is important to recall that 
the magnetism on a granular scale, which forms the basis for 
micromagentics, has its origins in atomic-level crystal field 
and exchange effects. The uniaxial anisotropy associated with 
hexagonal crystal systems gives rise to the well-known single- 
site contribution to the energy of the form -K{M*f. For K>0, 
this term is minimized by M||c||z, giving rise to the desired 
perpendicular anisotropy for crystallites having the c axis 
perpendicular to the film plane. However, hexagonal 
symmetry also allows for inter-site anisotropy terms of the 
form -J z (M i z Mj Z ). Such anisotropic exchange has been 
examined in the cases Co and Fe ions [2,3] and has previously 
been included in a micromagnetic calculation of domain 
structures in MnAs films [4]. 

The present work reports on an initial investigation of 
such anisotropic exchange in a micromagnetic model relevant 
to modern Co-based perpendicular recording media. 
Emphasis is placed on the impact of this interaction on the 
relative shapes of easy and hard-axis M-H loops. A specific 
application is made to the exchange coupled composite 
(ECC) media reported by Wang et al. [5]. The results of 
fitting the model to experimental M-H loops for both hard and 
soft layers individually are shown. Comparisons with data on 
loops in the case of exchange coupled films are also made. 
The effect of anisotropic exchange on thermal stability is 
discussed and some conclusions are also made regarding the 
benefits of ECC media. 



II. Micromagnetic simulations 

Simulations were performed using the commercial 
Landau-Lifshitz-Gilbert (LLG) micromagnetic simulator [6] 
with rectangular-prism grains on a square lattice with 32x32 
grains and periodic boundary conditions. Both domain-wall 
(incoherent) and rotational (coherent) reversal behavior can 
occur within this discretization scheme. Anisotropic exchange 
is incorporated as a modification of the usual isotropic 
exchange term by writing 

E ex = ~J Y.tj(MfMf + MjMj +mMfMf ) (1) 

where the sum is over near-neighbor grains, m>l induces 
uniaxial anisotropy and m<l creates planar anisotropy. Note 
that anisotropy of this type does not occur in crystals with 
cubic symmetry. For the exchange coupling between hard and 
soft layers in our model of ECC media, it was always assumed 
that m=l. The usual exchange parameter is expressed as 
A=Jd 2 , where d 2 is the cross-sectional area of a grain. Single- 
ion uniaxial anisotropy (K), fourth-order uniaxial anisotropy 
(K 2 ) [7], and magnetostatic interactions are also included. 
Zero temperature calculations were performed using a Suzuki- 
Trotter rotation-matrix method [8] with a damping parameter 
a=l. Finite temperature effects were included through the 
usual Langevin stochastic term [9] with an Euler integration 
routine, a=0.2, and a time step of 0.1 ps. Various quantities 
such as the saturation magnetization M s , anisotropy magnitude 
(K) and direction (k) were assumed to have a Gaussian 
distribution characterized by a given standard deviation, SD. 
We note that distributions in magnetic properties also serve to 
account for distributions in geometrical features, such as grain 
size (not explicitly included in the present model). 

III. IMPACT OF ANISOTROPIC EXCHANGE ON M-H LOOPS 

In order to explore the impact of anisotropic exchange on 
M-H loops, a generic model of single-layer perpendicular 
media was assumed. Grains with dimensions 7x7x20 nm 3 



BE-08 



2 



were used, having the following properties: M s =4.5xl0 5 A/m 
(SD=10%), ,4=0.05 - 0.15 xlO" 11 J/m, anisotropic exchange 
parameter m=l or 2, AM.OxlO 5 J/m 3 (SD=10%), as well as 
k sd =4°. Figs. 1 and 2 illustrate some key features, with the 
applied field H normalized by the anisotropy field H K =2K7M S . 
With m=2, both the coercivity and slope increase for out-of- 
plane (easy axis) loops compared to the isotropic case, m=l. 
Reducing the exchange parameter A has the effect to reduce 
the coercivity (H c ) and the slope, as expected [10]. In the case 
of in-plane loops, the slope decreases substantially for m>l 
with little effect on the coercivity. This is a noteworthy 
result. These features are consistent with the expectation that 
increasing m will serve to both increase the effective 
anisotropy as well as the exchange interaction if M z ^0. 
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Fig. 1. Easy-axis M-H loop illustrating some effects of 
anisotropic exchange in perpendicular media. The overall 
exchange parameter A is in units of 10" 11 J/m. 
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Fig. 2. As in Fig. 1 but with the applied field in the plane. 

IV. MODEL OF THE ECC MEDIA OF WANG ETAL. 

A. Hard Layer 

Fig. 3 shows the results of fitting model parameters for both 
cases with isotropic and very anisotropic (m=5), exchange in 
an attempt to reproduce the out-of-plane M-H loop data for the 
hard layer of Wang et al. [5]. Grains with dimensions 7x7x20 
nm 3 were again used along with M s =4.5xl0 5 A/m (SD=10%). 
For the case with m=5, fourth-order uniaxial anisotropy 
(K 2 =K/4) was also included with the same SD as assigned to 
K. A large easy-axis-direction SD of k sd =10° was required 
(for both m=l and m=5) to reproduce the rounding observed 
near the nucleation field (H n ). Larger k S d also has the effect of 
decreasing the easy-axis H c and increasing the hard-axis H c . 
This value of k sd required to reproduce the data is much larger 



than the reported A8 50 = 3° - 4°. The reasons for this difference 
are not entirely clear but may be related to grains having a 
core-shell structure, with a large distribution in the shell 
anisotropy axes [11]. In the case of m=l, the fitted values of 
exchange and anisotropy were found to be ^4=0.15xl0~ u J/m 
and A=2.0xl0 5 J/m 3 (SD=20%). For the m=5 case, the values 
^=0.02x10"" J/m and K=l. 45x10 s J/m 3 (SD=20%) gave the 
best fit to the data. 

Using these fitted parameters, somewhat better agreement 
with the in-plane experimental data is found in the case of 
m=5, as shown in Fig. 4. This can be understood since the 
exchange parameter A is reduced from its m=l fitted value, 
thus giving a smaller slope as seen experimentally. However, 
the overall agreement is significantly worse than in the easy- 
axis case of Fig. 3. Note the large experimental H c , not 
expected from coherent rotational switching usually assumed 
for hard-axis loops, as well as the lack of saturation in M. 
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Fig. 3. Model fit to the out-of-plane hard layer M-H loop data 
of Wang et al. [5] for both isotropic and anisotropic exchange. 
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Fig. 4. In-plane hard layer M-H loops resulting from using 
out-of-plane fitted parameters from Fig. 3. Experimental data 
from Wang et al. [5] is also shown. 

B. Soft Layer 

A reduction in film thickness as well as Pt content 
(relative to the hard layer) serves to increase the amount of 
stacking faults [5]. The effect is to reduce H c (by about a 
factor of 10) in the soft layer and produce a more isotropic 
film, with a larger anisotropy axis distribution. Model 
parameters were adjusted to best reproduce the experimental 
soft layer out-of-plane M-H loop [5]. Here, the layer 
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thickness was reduced to 8 nm and a saturation magnetization 
value M s =8xl0~ (SD=10%) was assigned. For the case m=l, 
the values k sd =17°, ^=0.48x10 _11 J/m, K=3.8xl0 5 J/m 3 
(SD=10%) fit the data well, as shown in Fig. 5. Assuming 
m=5, the fitting procedure yielded values k sd =20°, 
^=0.075x10"' 1 J/m, AT=1.25xl0 5 J/m 3 (SD=20%) which 
notably required larger distributions, which are known to 
degrade playback signal quality [12]. Note that the anisotropy 
fields, H K =2K/M S , are about the same for hard and soft films 
but the intrinsic anisotropy parameter D=K/M S 2 is smaller in 
the soft film by about a factor of two. The reduction in soft- 
layer He is mainly due to a larger M s and k S d- 

Using these parameter values, modeled in-plane loops, 
compared with the experimental data, are shown in Fig. 6. As 
with the hard layer, agreement with the data is much worse but 
again the model with anisotropic exchange appears somewhat 
better. Note that the experimental values of H c in both Figs. 5 
and 6 are similar, suggesting that a combination of domain- 
wall and rotational reversal occurs for both hard and easy -axis 
loops. The model again appears to underestimate the 
anisotropy magnitude (see Conclusions). 




-600 -400 -200 O 200 400 600 
H(kA/m) 

Fig. 5. Model fit to the out-of-plane soft layer M-H loop data 
of Wang et al. [5] for both isotropic and anisotropic exchange. 
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Fig. 6. In-plane soft layer M-H loops resulting from using out- 
of-plane fitted parameters from Fig. 3. 

C. Dual Layer 

ECC media is fabricated by controlled exchange 
coupling across soft and hard films through varying the 
thickness t s of a spacer layer, composed of Pt in the work of 
Ref. [5]. Parameters identified by the above analysis were 
used to model such coupled media. An inter-layer exchange 
parameter, A=Q.Q - 1.5xlO -11 J/m was added (with [13] m=l) 
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and varied in an attempt to mimic the effects of the 
experimentally controlled t s =0.7 - 3.6 nm. 

Some example easy-axis loops are shown in Figs. 7 and 8 
for the cases of m=l and m=5, respectively, where A t =0 and 
Aj = l.OxlO" 11 J/m were used. Data from Ref. [5] for the 
thinnest and thickest spacer layers are also shown. The 
general trends that H c is little affected, but that slopes decrease 
(increasing saturations fields H s ) with increasing t s , is 
reproduced by the model using a smaller A;. It is also evident 
that the model with anisotropic exchange again gives results in 
better agreement with the data. Overall, the observed effect of 
reduction in H c by about 40% from the single-layer value is 
reproduced by the model. Hard-axis loops were also 
calculated (not shown) which exhibit coercive fields of about 
70 kA/m for the m=l case and 25 kA/m for the m=5 case. 



Dual Layer: Out-of plane (m=l) 




Fig. 7. Easy-axis dual-layer media loops comparing m=l 
modeled and experimental results [5]. Interlayer exchange Aj 
is in units of 10" 11 J/m. 
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Fig. 8. As in Fig. 7 but with m=5 modeled results. 

V. Magnetization vs temperature. 

A useful indication of the impact of thermal fluctuations in 
magnetic systems is the reduction in M with increasing 
temperature [14]. Examination of simulation results for M vs 
time exhibit an equilibrium value. For the finite-size systems 
(as studied here), the transition to paramagnetism, T c , is 
indicated by a large reduction in M from its zero -temperature 
value (reaching zero only for infinite systems) [15]. 

Finite-temperature simulations were performed on the 
single and dual layers using model parameters discussed 
above to obtain the magnitude of M vs T. That thermal effects 
were properly treated in the model was verified by 
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reproducing the results of Ref. [14] as well as corresponding 
MC simulations. Simulation times of 150 ns were found to be 
adequate for the purpose of estimating M(T). We estimate 
(somewhat arbitrarily) T c by the relation M(T C ) = 0.15M(0). 

Results for the hard-layer magnetization, in single-layer 
cases m=land m=5, as well as the case of m=l for dual layers 
with Ai=0.1 and l.OxlO" 11 J/m, are shown in Fig. 9. In the 
single-layer cases, T c for both is estimated to be approximately 
1700 K. Interlayer exchange has little impact on T c when the 
hard layer is coupled to the soft layer. It is estimated to be 
about 1200 K for both dual-layer cases studied. A more 
accurate estimate of T c will require simulations using larger 
lattices. 
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Fig. 9. Hard layer magnetization as a function of 
temperature for single and dual layer films. Interlayer 
exchange A t is in units of 10" 11 J/m. 

VI. SPIN WAVES 

If magnetic interactions in sputtered cobalt -based films do 
include significant exchange anisotropy, it is not immediately 
obvious how this can be measured. One possibility is through 
spin-wave dispersion curves as a function of wave vector q. 
Single-ion anisotropy K appears as a gap at q=0 and the curve 
shape determines the exchange constant J. A simple 
calculation based on the torque equation (LLG equation with 
a=0) shows that for perpendicularly magnetized thin films, K 
and f simply add together to increase the gap [16]. However, 
if a strong magnetic field H is applied in the film plane, 
anisotropic exchange also affects the shape of the dispersion 
as in the following relations (for a square lattice with inter- 
atomic separation distance a): 

a> 2 q = [H + 4Af/(l + f q )][H - 2KM + 4M/(l - mf q )] (2) 

f q — (cosq x a + cosq y a)/2 (3) 

The anisotropy parameter m could in principle be used as a 
distinguishable fitting parameter along with J and K. 

VII. DISCUSSION AND CONCLUSION 
This work has shown that anisotropic exchange may play a 



role in micromagnetic models of perpendicular media, 
especially when faced with the challenging task of 
reproducing both easy and hard axis M-H loops. We note that 
an approach based on mimicking long experimental sweep- 
rate effects by assuming an elevated temperature [17] yields 
larger fitted anisotropy values and better agreement for in- 
plane loops (to be reported elsewhere). 

Some conclusions from the present work may be made 
regarding the relative benefit of ECC media. Although there 
is a reduction in H c , we have found a substantial reduction in 
T c as well when the hard layer is coupled to a soft layer 
(indicating larger thermal fluctuations). Taking the ratio of 
H c /T c as a figure of merit, our results yield approximately 3.0 
kA/(mK) for the single-layer and 2.5 for the dual-layer cases, 
where a smaller number is advantageous for a combination of 
writability and thermal stability. The modeled ECC media is 
seen to yield about a 20% advantage by this criterion. 
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